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The regulation of plant cell wall synthesis pathways remains poorly understood. This has become a 
bottleneck in designing bioenergy crops. The goal of this study was to analyze the regulation of plant cell 
wall precursor metabolism using metabolic flux analysis based on dynamic labeling experiments. 
Arabidopsis T87 cells were cultured heterotrophically with 13 C labeled sucrose. The time course of 13 C 
labeling patterns in cell wall precursors and related sugar phosphates was monitored using liquid 
chromatography tandem mass spectrometry until steady state labeling was reached. A kinetic model 
based on mass action reaction mechanisms was developed to simulate the carbon flow in the cell wall 
synthesis network. The kinetic parameters of the model were determined by fitting the model to the 
labeling time course data, cell wall composition, and synthesis rates. A metabolic control analysis was 
performed to predict metabolic regulations that may improve plant biomass composition for biofuel 
production. Our results describe the routes and rates of carbon flow from sucrose to cell wall precursors. 
We found that sucrose invertase is responsible for the entry of sucrose into metabolism and UDP- 
glucose-4-epimerase plays a dominant role in UDP-Gal synthesis in heterotrophic Aradidopsis cells 
under aerobic conditions. We also predicted reactions that exert strong regulatory influence over carbon 
flow to cell wall synthesis and its composition. 

© 2013 Elsevier Inc. All rights reserved. 


1. Introduction 

Plant cell wall material is a major component of plant biomass 
that mostly consists of sugar polymers (Somerville et al., 2004). It 
is the most abundant sugar resource on earth and a promising 
feedstock for biofuel production (Pauly and Keegstra, 2010). One 
major bottleneck in plant biomass based biofuel production is the 
resistance of plant cell wall to biological, chemical or physical 
deconstruction (Chundawat et al., 2011). Understanding plant 
synthesis pathways and their regulation may lead to break¬ 
throughs in bioenergy crop design. Current research on plant cell 
wall synthesis is mainly focused on understanding how cell wall 
precursors were assembled to polysaccharides at the gene and 
protein levels (Bar-Peled and O’Neill, 2011; Carpita, 2012; Keegstra 
and Cavalier, 2011; Martinez et al., 2009; Zhang et al., 2011). The 
literature indicates that plant cell wall synthesis has not yet been 
predictively modeled at the metabolic network level. Metabolic 
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flux analysis (MFA) provides the necessary tools to study the 
patterns of carbon flow and their regulation, which underpin plant 
growth and development. 

In MFA internal carbon fluxes through metabolic networks are 
obtained using isotope labeling experiments (Antoniewicz et al., 
2007; Nanchen et al., 2007; Sauer et al., 1999; Schmidt et al., 1997; 
Wiechert and de Graaf, 1997; Wiechert et al., 1997; Wittmann, 
2002; Zamboni et al., 2005, 2009). Most MFA studies in plant 
research are performed at both isotopic and metabolic steady 
state, where labeled substrates are supplied to plant tissues/cells, 
carbon rearrangements of metabolites are tracked and fluxes are 
estimated by optimized fitting of internal fluxes to the measure¬ 
ments of labeling and external fluxes (substrate uptake and 
biomass production rates). It has been successfully applied in 
different plant systems to quantify carbon flow through central 
metabolism (for recent reviews see Allen et al., 2009; Chen and 
Shachar-Hill, 2012; Kruger et al., 2012; Kruger and Ratcliffe, 2009; 
Libourel and Shachar-Hill, 2008; Schwender, 2011). However, 
steady state MFA is limited to metabolic steps and pathways that 
involve rearrangement of the label (usually carbon). These limita¬ 
tions lead to the need for dynamic metabolic flux analysis. 

Dynamic metabolic flux analysis (DMFA) is an MFA method 
based on transient labeling experiments (for reviews see Noack 
et al., 2011; Noh and Wiechert, 2011; Wiechert et al., 2005). It 
involves tracking the propagation of labeling patterns through a 
metabolic network and fitting a kinetic model (in which ordinary 
differential equations are used to represent reaction rates) to the 
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dynamic labeling measurements obtained from pulse labeling 
experiments (without chase or washout). It is able to describe 
metabolic activities of a network in vivo and predict the metabolic 
effects on network fluxes in response to perturbations. Despite its 
advantages over steady state MFA, DMFA often requires the 
development of case specific analytical methods for metabolic 
intermediate measurements as well as the use of mathematical 
modeling. Computational frameworks for DMFA are emerging and 
have been demonstrated or applied in different biological systems 
(Leighty and Antoniewicz, 2011; Noh et al., 2006; Selivanov et al., 
2004, 2006; Young et al., 2008), but none of these frameworks has 
been implemented as widely applicable software. Due to these 
difficulties, DMFA in plant research has been mostly performed in 
small to medium sized metabolic networks that have little or no 
carbon skeleton rearrangements (Boatright et al., 2004; Colon et al., 
2010; Knoke et al., 2009; Matsuda et al., 2003, 2005, 2007; McNeil 
et al., 2000; Rios-Estepa et al., 2008; Rohwer and Botha, 2001). 

Here we have developed a DMFA model to study the regulation 
of plant cell wall precursor metabolism. We have monitored the 
time course of 13 C labeling in cell wall precursors and related sugar 
phosphates using a liquid chromatography tandem mass spectro¬ 
metry method that we previously developed (Alonso et al., 2010). 
We used the DMFA model based on mass action reaction mechan¬ 
isms to simulate the carbon flow in the cell wall synthesis 
network. The kinetic parameters of the model were determined 
using the time course labeling data of the cell wall precursors. The 
model allowed us to explore key reactions that control carbon flow 
to cell wall synthesis and predict in vivo metabolic regulation that 
may affect plant cell wall composition. 


2. Materials and method 

2.2. Chemicals 

UDP-Xyl, UDP-Ara and UDP-GalA were purchased from the 
Complex Carbohydrate Research Center (University of Georgia). 


All the other NDP-sugars and hexose phosphates, excluding UDP-Api 
and UDP-Rha which are not commercially available, were purchased 
from Sigma. [U- 13 C Fm ]Suc and [U- 13 C G i c ] were, respectively, pur¬ 
chased from Omicron Biochemicals, Inc and Isotec (Sigma-Aldrich). 

2.2. Arabidopsis cell culture, sampling and extraction of metabolites 

Arabidopsis thaliana ecotype Columbia (cell line T87,Jouanneau 
and Peaud-Lenoel, 1967) were obtained from the Riken BioRe- 
source Center. Cell culture and sampling were performed as 
previously described (Alonso et al., 2010). Briefly, Arabidopsis 
T87 cells were cultured in a 7 day culture interval under contin¬ 
uous dark conditions at 22 °C and 120 rpm in a liquid minimal 
media (Alonso et al., 2010). T87 cells at log phase were centrifuged 
at 2000 rpm for 3 min at 22 °C. Supernatant was removed under 
sterile conditions. Cells were transferred to the fresh minimal 
media containing 100% 13 C labeled sucrose ([U- 13 C g i c ]-sucrose or 
[U- 13 C fm ]-sucrose). Cells from 1 mL suspension culture were 
sampled by fast filtration via a nylon membrane filter (Whatman, 
0.45 pm pore size and 47 mm diameter) and frozen in liquid 
nitrogen. Intracellular metabolites were extracted using boiling 
water (Alonso et al., 2010). The time course of 13 C labeling in 
cell wall precursors and related sugar phosphates were monitored 
using a liquid chromatography tandem mass spectrometry method 
(Alonso et al., 2010) until isotopic steady state was reached 
(Supplementary materials Tables S3 and S4). To measure cellular 
sucrose, fructose and glucose concentration, 3 mL cell suspension 
was filtered through the membrane filter and the filtrates on the 
membrane were washed with 6 mL media lacking sucrose. Cellular 
sugars were measured using Megazyme® sucrose/fructose/ 
D-glucose kits. 

2.3. External flux measurements 

T87 cell growth was measured by monitoring the OD 60 o nm- Cell 
suspension was collected at different time points. Cells and cultured 
medium were separated by centrifugation at 13,000 rpm for 5 min 


*Gal1P 



Fig. 1 . Metabolic network model of plant cell wall synthesis. Bold font (dash line) represents concentration (flux) measurements were made. Star represents labeling 
measurements were made. Italic font represents concentration measurements that include indistinguishable quantities (GDP-Man and GDP-Glc; and Ml P and G1P) and were 
measured as GDP-HEX and HIP, Dot line represents flux measurements that include indistinguishable quantities (GDP-Glc -> Glc and UDP-Glc -> Glc). Total glucose 
monomer content in cell wall was used to obtain the sum of two effluxes that depose glucose into cell wall. Circle represents the measured labeling was treated as inputs to 
the model. Boxes represents terminal of carbon flow. 
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at 4 °C. Cells were washed with deionized water twice and 
lyophilized. Glycosyl composition analysis of cell wall was per¬ 
formed at the Complex Carbohydrate Research Center by combined 
gas chromatography/mass spectrometry of the per-O-trimethylsilyl 
derivatives. Starch measurements were performed as previously 
described (Alonso et al., 2007), with lipid extraction, followed by 
washing with 1 mL 80% ethanol twice and starch extraction, 
hydrolysis and enzyme-based glucose measurement. The starch 
concentration in cell debris was measured using the Megazyme® 
Starch Kit. Extracellular uptake rates and production secretion 
rates were calculated using the equation (Sauer et al., 1999): 
v = ((C(t 2 )-C(ti))/(OD(t 2 )-OD(ti))-/) y M, where C(t) and OD(t) are 
the extra-cellular metabolite concentrations and OD 60 o at time t; 

is the growth rate of the culture; and/is the conversion coefficient 
between OD 60 q and biomass dry weight per unit volume of media. 


Table 1 

Growth rate, external flux and intracellular metabolite measurement. 


f* (1/day) 

0.41 

± 

0.04 

/ (g dry weight/od/1) 

2.97 

± 

0.62 

Cell wall/ Cell biomass dry weight 

18% 

+ 

1% 

Csucrose (nmol/mg dry weight) 

232.4 

± 

30.3 

Cgiucose (nmol/mg dry weight) 

50.5 

± 

10.1 

Cfructose (nmol/mg dry weight) 

9.1 

± 

5.1 

V S uc_ex (nmol/mg dry weight/day) 

2561.4 

+ 

102.3 

V G ic_ex (nmol/mg dry weight/day) 

445.3 

+ 

168.7 

Vstarch (nmol/mg dry weight/day) 

1.0 

± 

0.4 

V Rha (nmol/mg dry weight/day) 

21.8 

± 

2.0 

V Fuc (nmol/mg dry weight/day) 

4.1 

± 

1.3 

V Ara (nmol/mg dry weight/day) 

57.3 

± 

3.9 

V Xy i (nmol/mg dry weight/day) 

21.9 

± 

0.5 

Vivian (nmol/mg dry weight/day) 

9.1 

± 

1.1 

V Ga i (nmol/mg dry weight/day) 

49.8 

± 

2.1 

V G ic (nmol/mg dry weight/day) 

208.1 

± 

12.3 

V Rib (nmol/mg dry weight/day) 

10.5 

± 

0.3 

Vgica (nmol/mg dry weight/day) 

0.7 

± 

0.3 

VcaiA (nmol/mg dry weight/day) 

71.7 

± 

14.0 

Vcic (nmol/mg dry weight/day) 

208.1 

± 

12.3 


2.4. Metabolic network model of plant cell wall precursors 


A metabolic network of plant cell wall precursors were con¬ 
structed based on previous pathway studies (Alonso et al., 2010; 
Reiter, 2008; Seifert, 2004; Sharpies and Fry, 2007) (Fig. 1). We 
developed a DMFA model based on the mass balance principle to 
simulate the flow of carbon cumomers in the network: 
d(m[-Ci)/dt = ZR J m \ where Q is the concentration of the precursor 

i, m[ is the cumomer m k of precursor i, and Rj^is the flux of the 

reaction j that consumes or produces precursor i with cumomer m k 
(/<=0,3,6). Here mass action kinetics were used to simulate fluxes 
through each reaction in the network. For example, the balance of 
carbon cumomer for GallP during dynamic labeling experiments 
(pulse labeling without a chase or washout) can be written as: 


^(^GallP’^Gallp) 

dt 


-^GallP UGal -^Gallp-^GallP G1P ■ m GallP 


+^GlP Gall p-^Glp 



The rate equations were listed in Table 1. 

The DMFA model was simulated using the PottersWheel soft¬ 
ware package (Maiwald and Timmer, 2008). The model was 
simulated in the form of customized ordinary differential equa¬ 
tions in the software. In addition to the equation system, two types 
of constraints were used in the simulation: (1) no change in the 
pool size of internal metabolites during the dynamic labeling 
experiments dCf/dt = ZRj = 0. This ensures the mass balance of 
the model in each time step of the parameter fitting; and (2) the 
ratio of carbon flow into pentose phosphate pathway and glyco¬ 
lysis pathway was constrained to be greater than 39% and less 
than 65% based on the previous 13 C metabolic flux analysis of 
heterotrophic Arabidopsis cell central metabolism (Masakapalli 
et al., 2010; Williams et al., 2008). 

Simulation were performed in a Dell Studio XPS 9100 PC 
(8 Intel l! Core™ CPU at 2.8 GHz, 8.0 GB RAM and 64-bit operating 
system). Starting from a set of initial parameters, the PottersWheel 
software gradually adapted the parameter in the ODEs to fit the 



Fig. 2. Intracellular metabolite levels. The concentrations of phosphate sugars and nucleotide sugars were obtained from Alonso et al. (2010). Cellular sucrose, glucose and 
fructose measurements were performed in this work. 
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labeling, concentration and flux measurements (minimizing the 
sum of squared residuals weighted by the standard deviation of 
each measurement (SS res )). 320 sets of initial parameter values 
were used here to avoid finding only local minima. The fittings 
were performed with disturbed initial parameter values using 
Trust Region Algorithm in log space. The best fitting result was 
used to provide the optimized rate constants. 

The confidence intervals of the optimized rate constants were 
estimated using a Monte Carlo approach (Antoniewicz et al., 2006; 
Libourel and Shachar-Hill, 2008). Fifty datasets were generated using 
normal distributions N(p } o) for each measured quantity, where g is the 
mean of the labeling, concentration or flux measurement and o is the 
standard deviation of the measurement. The generated data form 50 
sets of Monte Carlo pseudo-measurements representing random 
variations of the original measurements. The DMFA model was fitted 
to these pseudo-datasets as described in the previous section (104 
initial parameter sets were used for each set). The confidence limits 
were given by the standard deviation of the best fit optimized 
parameters for each of the 50 pseudo-datasets. 

Metabolic control analysis was performed using the software 
package COPASI 4.7 (Build 34) (Mendes et al., 2009). 


3. Results and discussion 

3.2. Model development 

Fig. 1 shows the literature derived network (Alonso et al., 2010; 
Reiter, 2008; Seifert, 2004; Sharpies and Fry, 2007) through which 
carbon flow from the substrate sucrose to cell wall precursors, the 
oxidative pentose phosphate pathway (OPPP) and glycolysis, and 
the flow from precursors into cell wall components. The formation 
of cell wall polymers was assumed to be terminal in the current 
model. Thus the assembly, modifications and any turnover of cell 
wall components after monomer polymerization were not con¬ 
sidered here. The OPPP and glycolysis pathway were treated as net 
sinks for carbon flow. 

The entry of label into both central metabolism and cell wall 
precursor pools, originated from the feeding 13 C labeled sucrose. 
However since there are reversible reactions in central metabo¬ 
lism, the cell wall network also received 13 C labeling from central 
metabolism via the common metabolite pools (glucose-6-phos- 
phate (G6P) and fructose-6-phosphate (F6P)) shared by central 
metabolism and the cell wall synthesis network. We therefore 



Fig. 3. Cell wall glycosyl composition (grams glycosyl residue/100 g biomass dry 
weight). 


used sucrose, G6P and F6P as labeling input points to the cell wall 
network in our modeling. 

Fig. 1 also summarizes the concentration, labeling and fluxes 
measurements that we conducted. We measured the concentration 
of 19 metabolites in the network and time course 13 C labeling of 10 
metabolites (M0, M3 and M6) in the network (Supplementary 
materials). We were not able to separately measure levels or 
labeling in two isomeric pairs of metabolites (ManlP and GlclP; 
and GDP-Man and GDP-Glc). The two metabolites were modeled 
separately, and the predicted average labeling was used to match 
the measurements during fitting. 

3.2. Fitting and simulation results 

320 sets of initial parameter values (rate constants and metabo¬ 
lites) were used here to obtain the global best fitting. The distribution 
of SS res from 320 fittings was shown in Fig. 4B. The global SS res /N is 
64.3 m, where N is the number of measurements used in the fitting 
(N=361). The optimized parameter values and their confidence 
interval in the form of standard deviations are shown in Table 2 and 
Supplemental materials Table SI. The mean of best fit parameters 
from Monte Carlo simulation were shown in Supplemental materials 
Tables SI and S2. Overall, 56% of the normalized (to the mean of 
optimized kinetic constants) standard deviation is smaller than 0.3. For 
the normalized standard deviations that were greater than 30%, 35% of 




Fig. 4. (A) Correlation plot between measurement and prediction; and (B) fitting 
residues from 320 fittings. 
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Table 2 

Rate equation and optimized rate constant. 


Flux Name 


1 

V uptake 

2 

V S uc synthase 
V S uc synthase 

3 

V S uc invertase 

4 

Vf6P-M6P 

Vm6P-F6P 

5 

Vm6P-M1P 

Vm1P-M6P 

6 

VMIP-GMan 

WMan-MlP 

7 

Vefflux-Man 

8 

V efflux-Fuc 

9 

VGallP-UGal 

10 

VGallP-GlP 

WlP-GallP 

11 

VuGal-UGlc 

VuGlc-UGal 

12 

WlP-UGlc 

13 

Vg1P-G6P 

Vg6P-G1P 

14 

VuGlc-URha 

15 

VuGlc-UGlcA 

16 

VuGlcA-UGalA 

VuGalA-UGlcA 

17 

Vg6P-6PG 

18 

Vefflux-Gal 

19 

V e fflux-Glcl 

20 

Vgip-ggic 

Vggic-gip 

21 

V efflux_Glc2 

22 

V e fflux_Rha 

23 

V efflux_xyl 

24 

V efflux_Ara 

25 

V e fflux_GalA 

26 

V e fflux_OPPP 

27 

V efflux_ex_Glc 

28 

V efflux_starch 

29 

Vfi-u-F6P 

30 

Vg1c-G6P 


Direction 


Reaction 


Suc ex — > SUCcellular 
SucoUDP-Glc+Fru 

Suc->Glc+Fru 

F 6 P 0 M 6 P 

M 6 P 0 MIP 

MIP^GDP-Man 


GDP-Man^ 

GDP-Man^ 

Gall UDP-Gal 
GallPoGlP 


UDP-GaloUDP-Glc 

GIP^ UDP-Glc 
GIP 0 G 6 P 

UDP-Glc -► UDP-Rha 
UDP-Glc ->■ UDP-GlcA 
UDP-GlcAoUDP-GalA 


G 6 P-> 6 PG 

UDP-Gal^ 

UDP-Glc^ 

GlPoGDP-Glc 

GDP-Glc^ 

UDP-Rha^ 

UDP-GlcA^ 

UDP-GlcA^ 

UDP-GalA^ 

6PG^ 

Glc^ 

GIP^ 

Fru -> F 6 P 
Glc->G 6 P 


Rate equation 

Rate constant (k) 

std 

K 

2.56E+03 

9.85E+1 

K*C Suc 

1.80E-02 

5.57E-03 

K*CuDP-Glc*Qru 

1.52E-05 

2.04E-03 

I<*C suc 

1.13E+01 

4.38E-01 

ICC F6 P 

8.93E+00 

7.42E-01 

ICC M6 p 

9.65E-01 

5.43E-01 

ICC M 6p 

1.83E+01 

1.48E+00 

I<*Cmip 

5.01 E-01 

1.57E+00 

K*C M ip 

7.48E+01 

1.76E+01 

K*C G DP-Man 

1.30E+01 

7.77E+01 

K CcDP-Man 

3.33E+02 

4.37E+01 

K*C G DP-Man 

1.26E+02 

2.16E+01 

I<*C G allP 

6.97E-01 

1.31E+00 

I<*C G allP 

5.33E-02 

2.18E-01 

ICCgip 

5.37E-01 

9.42E-01 

K*CjjDP-Gal 

1.69E+00 

1.61 E+00 

K*CuDP-Glc 

3.41 E+01 

1.41 E+00 

K*C G ip 

2.97E+03 

9.67E+02 

ICCgip 

6.17E+00 

1.00E+01 

ICCgsp 

6.36E+01 

3.03E+00 

ICCudP-GIc 

1.11 E+01 

1.73E+00 

iCCuDP-Glc 

9.82E+01 

9.32E+00 

ICCudP-GIcA 

6.77E+02 

1.23E+02 

iCCuDP-GalA 

3.10E+01 

1.18E+01 

ICCgsp 

2.68E+02 

2.70E+01 

ICC UDP-Gal 

1.47E+02 

5.76E+00 

ICCuDP-Glc 

1.40E+02 

8.44E+00 

ICCgip 

3.44E+00 

3.38E+00 

iCCcDP-Gld 

5.98E+01 

3.90E+01 

ICCgdp-gic 

4.51 E-01 

1.16E+01 

iCCuDP-Rha 

1.64E+01 

2.56E+00 

ICCudP-GIcA 

2.22E+02 

1.81 E+01 

ICCudP-GIcA 

4.63E+02 

5.54E+01 

ICCuDP-GalA 

2.68E+03 

4.96E+02 

ICCspg 

1.57E+06 

1.58E+05 

ICCgic 

4.02E+00 

1.74E+00 

ICCgip 

7.42E+00 

2.92E+00 

iCCpru 

3.07E+02 

1.18E+01 

ICCgic 

4.91 E+01 

2.77E+00 


them came from reactions with small fluxes. A correlation plot of the 
measurement constraints used in fitting and model predicted results is 
shown in Fig. 4A. The labeling prediction for several metabolites 
during the dynamic labeling experiments was shown in Fig. 5. 

3.3. Physiology 

Figs. 2 and 3, and Table 1 depict the physiology of T87 cells in the 
culture conditions used. The cells were maintained in metabolic steady 
state (the pool sizes of metabolic intermediates remained the same 
through the time course) during the dynamic labeling experiments. 
T87 cells grew in the current culture condition at a growth rate of 
0.41 ± 0.04 day -1 . The glycosyl components in the cell wall accounted 
for 18% +1% of total biomass dry weight. This corresponds to a 
significant portion of carbon uptake being incorporated into the cell 
walls (-8.8%). This is in the range of previous findings in Arabidopsis 
cell culture (Masakapalli et al., 2010; Williams et al., 2008). 

3.4. Flux map 

The flux map of cell wall synthesis network at metabolic steady 
state is shown in Fig. 6. The results describe the routes and rates of 
carbon flow from sucrose to cell wall precursors. We found that 
sucrose invertase is responsible for the entry of sucrose into 
metabolism in T87 cells. More than 99% of carbon entered cellular 
metabolism via sucrose invertase under these growth conditions. 
The carbon flow was then directed to cell wall precursor synthesis 
through phosphoglucomutase and phosphomannose isomerase. 


Our results provide evidence on the role of sucrose invertase in 
cellulose and starch production, which has been a topic of debate 
in recent years (Baroja-Fernandez et al., 2012; Barratt et al., 2009; 
Bieniawska et al., 2007). Barratt et al. (2009) generated Arabidop¬ 
sis mutant plants that lack isoforms of sucrose invertase and 
sucrose synthase, respectively. It was found that the plants that 
lack sucrose synthase activities had normal growth phenotypes, 
but the sucrose invertase mutant plants had severely reduced 
growth rates. Baroha-Fernandez et al. used a different method to 
assay sucrose synthase activity in plant extracts and reported that 
the activities can support normal starch and cellulose biosynthesis. 
Their work, however, did not assay the activities of sucrose 
invertase in the same plant extracts. Our results provided evidence 
that in vivo invertase activity carries about over 99% of the entry of 
sucrose into metabolism. This supports an important role for 
sucrose invertase in the entry of sucrose into metabolism in 
heterotrophic cells under aerobic conditions. Our results suggested 
that phosphoglucomutase (Glc6P->GlclP) may act as a critical 
step that directs carbon flow to cell wall synthesis. 

Our flux map shows that UDP-Glc-4-epimerase plays a domi¬ 
nant role in UDP-Gal synthesis in heterotrophic Arabidopsis cells 
culture. UDP-Gal can be synthesized via unidirectional UDP-Gal 
pyrophosphorylase and bidirectional UDP-Glc-4-epimerase in 
plants. Sharpies and Fry conducted a dual-radiolabeling experi¬ 
ment on Arabidopsis cells cultured in a glycerol medium (Sharpies 
and Fry, 2007). The cell culture was fed with traces of D-[1- 3 H] 
galactose and [U- 14 C] fructose. 3 H: 14 C ratios decreased in the order 
of GallP, UDP-Gal, UDP-Glc, G1P, G6P and F6P. These results 
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Fig. 5. Labeling measurements and model prediction. (A) GallP in [U- 13 C fm ]-sucrose experiments; (B) GallP in [U- 13 C fm ]-sucrose experiments; (C) HexlP in [U- 13 C fru ]- 
sucrose experiments; and D: HexlP in [U- 13 C fm ]-sucrose experiments. 


GallP 


UDP-Gal 


Sue 



Fig. 6. Flux map at metabolic steady state. Flux unit: nmol/mg DW. 


demonstrated the existence of both UDP-Gal pyrophosphorylase 
and UDP-Glc-4-epimerase activities in the Arabidopsis cell culture. 
The relative importance of the two enzymes in UDP-Gal synthesis 
was not determined in that work. Our results show that UDP-Gal is 
almost solely synthesized from UDP-Glc via UDP-Glc-4-epimerase. 
This finding is also consistent with previous reverse genetic 
studies. Rosti et al. (2007) studied A. thaliana mutant plants and 


found that UDP-Glc-4-epimerase isoforms cooperate in providing 
UDP-Gal for cell wall biosynthesis and growth of A. thaliana . The 
dominant role of UDP-Glc-4-epimerase over UDP-Gal pyropho¬ 
sphorylase suggests that the former enzyme may be a more 
promising target for designing biofuel crops. 

We considered whether inaccuracies in pool sizes might affect 
the flux map and found that the fluxes derived from the modeling 
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Table 3 

Reactions that may affect carbon flow to cell wall xylose. 


Reaction 

Flux control coefficient^) 

Glc6P^GlclP 

99 

UDP-Glc -> UDP-GlcA 

65 

UDP-GlcA -> UDP-GalA 

-49 

UDP-Glc^cell wall Glc 

-50 

UDP-Ara -> cell wall Ara 

-34 


are independent of the pool sizes used. This can be shown both 
theoretically and empirically. Theoretically, each flux of material 
(reaction rate) is modeled as the product of a substrate concentra¬ 
tion and a rate constant, and each flux of label (rate of 13C flow) is 
modeled as the product of the product of the same substrate 
concentration and rate constant multiplied by the labeling level of 
the substrate. Therefore the value of a metabolite concentration 
always appears multiplied by the relevant rate constant. Therefore 
if a metabolite concentration was incorrect because of compart- 
mentation, modeling with a corrected concentration would result 
in the value of the rate constant obtained from the model being 
changed by the inverse of the same factor, leaving the flux 
unchanged. Empirically, we confirmed this by arbitrarily changing 
several of the concentrations used in the modeling and observing 
that the change of the flux values obtained are all within the 
confidence intervals calculated from the Monte Carlo approach 
(less than 1.9% of total carbon uptake rates) although the rate 
constants were different as predicted. 

3.5. Predictions 

One important goal of this study was to understand the 
regulation of cell wall synthesis network aid the design of plants 
with improved cell wall compositions. Here we performed meta¬ 
bolic control analysis to obtain insights into the regulation of 
cell wall precursor synthesis (Burns et al., 1985; Fell, 1992, 1997). 
The resulting flux control coefficients (Supplementary materials 
Table S5) depict how modest perturbations of enzyme activities 
can affect metabolic fluxes in the cell wall synthesis network. 

Pentose sugar contents such as xylose in plant cell wall are 
difficult to ferment by microbes, therefore representing an impor¬ 
tant target for metabolic engineering of bioenergy crops. Table 3 
highlighted the perturbations that may significantly affect carbon 
flow to cell wall xylose synthesis. The results show that doubling 
the expression of phosphoglucomutase and UDP-Glc dehydrogen¬ 
ase (UDP-Glc ^UDP-GlcA) may increase the carbon flux to cell wall 
xylose by 99% and 65%, respectively, and doubling the expression 
of UDP-GlcA-4-epimerase (UDP-GlcA ^UDP-GalA) and enzymes 
for cell wall glucose (UDP-Glc^ cell wall Glc) and cell wall 
arabinose (UDP-Ara->cell wall Ara) synthesis may reduce cell wall 
xylose formation significantly. Interestingly, the perturbations of 
enzymes that are adjacent to carbon entry (sucrose invertase) have 
stronger impacts on the cell wall xylose formation fluxes. 

3.6. Limitations 

The mathematical model developed here is presented as a tool 
for exploring the kinetic regulation of plant cell wall synthesis 
in vivo. The model was able to explain the time course data for 
multiple isotopically labeled forms (mass isotopomers) of multiple 
metabolites in terms of pathway delineation and prediction of 
metabolic regulation under steady state growth. However signifi¬ 
cant simplifications were made in the current approach. 

The network was modeled as a one-compartment system while 
some metabolites such as G6P and F6P appear in both cytosol and 


plastid. This was done because we were not able to reliably 
quantify levels or labeling of the same metabolite in different 
subcellular compartments. Steady state MFA studies have used 
compartment-specific products (Allen et al., 2007) to measure the 
labeling of sugar phosphate and nucleotide sugars in the cytosol 
and chloroplast (Lonien and Schwender, 2009; Masakapalli et al., 
2010; Williams et al., 2008). However steady state labeling 
patterns cannot reliably predict potential compartmental differ¬ 
ences during dynamic labeling. To overcome this challenge and 
refine the model it may be feasible to use non-aqueous fractiona¬ 
tion (Farre et al., 2001; Krueger et al., 2011) to isolate and analyze 
pools from cytosol and plastid. 

The kinetics were modeled as linear mass action rates for each 
reaction in the network. This simplification is because the dynamic 
labeling experiments were performed at metabolic steady state, 
therefore the labeling data cannot determine the non-linear 
enzymatic behavior or catalytic effects due to inhibitor/activator/ 
allosteric regulators. This approach is usual among dynamic 
labeling studies of networks of moderate or larger size (Libourel 
and Shachar-Hill, 2008). 

To relieve these limitations, methods for quantifying metabo¬ 
lites and cosubstrates at the subcellular compartment level need 
to be further developed (especially for labile metabolic intermedi¬ 
ates), and metabolomic profiling experiments at metabolic non¬ 
steady state will need to be conducted. 

4. Conclusion 

We have developed a kinetic model based on mass action 
reaction mechanisms to simulate the carbon flow in the cell wall 
synthesis network. We have determined the kinetic parameters of 
the model using the time course labeling data of the intermediates 
within the metabolic network. The results revealed sucrose 
invertase as the primary entering point of carbon into metabolism 
under non-photosynthetic cell suspension culture condition, and 
suggested that phosphoglucomutase (Glc6P^GlclP) may have a 
significant regulatory role in carbon flow to cell wall synthesis. 
The computational modeling also predicted in vivo metabolic 
regulation that may affect plant cell wall composition. 
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